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The interest in supersonic retro-propulsion (SRP) as a means of deceleration during 
planetary entry increases with the desire to land high mass vehicles on Mars. Since it is 
difficult to obtain flight data or properly simulate this type of flow field in a wind tunnel, the 
use of computational fluid dynamics (CFD) becomes increasingly important, as does the 
need to verify the current CFD methods. This presentation will show results from structured 
overset grids and OVERFLOW, a Reynolds Averaged Navier-Stokes solver, obtained during 
the continuing CFD verification process. Flow structure, surface pressure, forces, and 
moments are compared to historic and modern wind tunnel data as well as to other Navier- 
Stokes solvers, DPLR and FUN3D. Cases include single and multiple nozzle cases from the 
Jarvinen and Adams experiment, 1 the Daso et al experiment, 11 and a recent test in the NASA 
Langley Unitary Wind Tunnel (scheduled for June 2010). 


I Jarvinen, P.O. and Adams, R.H., “The Aerodynamic Characteristics of Large Angled Cones with Retrorockets,” 
NASA Contract No. NAS7-576, Cambridge, MA, Feb. 1 970. 

II Daso, E. et al, “Dynamics of Shock Dispersion and Interactions in Supersonic Freestreams with Counterflowing 
Jets,” AIAA Journal Vol. 47, No. 6, June 2009. 
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Introduction 


• The goal is to softly land high mass vehicles ( I Os of metric tons) on Mars 

• Supersonic Retro-Propulsion (SRP) is a viable means of deceleration 

• CFD is of increasing importance since flight and experimental data at 
these conditions is difficult to obtain 

• CFD must be validated at these conditions 






CFD Validation Approach 


• Employ multiple CFD codes (OVERFLOW, FUN3D, and DPLR) to solve 
the same SRP problems 

• New additions of US3D and Cart3D 

• Compare results between codes and with historic data 

• Shock structure/ shock standoff distance 

• Surface pressures and forces 

• Run new wind tunnel tests for CFD validation 

• Complete run conditions 

• Higher thrust coefficients than in existing data 
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Comparison 


Validation 
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Presentation Preface 


• Focus on OVERFLOW results 

• Results from Kerry Trumble (NASA 
ARC, DPLR) and Bil Kleb (NASA 
LaRC, FUN3D) will also be shown 

• This is a work in progress 
Outline 

1. Jarvinen and Adams Single Nozzle 

2. Daso et al Single Nozzle 

3. Jarvinen and Adams Triple Nozzle 

4. Langley UPWT 4x4 Pre-Test 

* Single, triple, and quad nozzle cases 

5. Current Work 

6. Future Work 


Langley UPWT pre-test, 4 nozzle, Mach 4.6, Cj 5 


Single Jet SRP Flow Field Structure 
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Freestream 
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Single Jet SRP Flow Field Structure 




Single Jet SRP Flow Field Structure 
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Single Jet SRP Flow Field Structure 
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Jarvinen and Adams Single Nozzle 




60 degree sphere cone 

Ames 6’x6’ supersonic wind tunnel 
in 1970 

Inconsistencies in the report 

• Geometric dimensions 

• Freestream total temperature and 
pressure 

• No uncertainties reported 

For this test, code-to-code 
comparison was relied on heavily 

Run conditions 

• Mach = 2.0, Re/in = 40604.3 

• C T = 0, 7 


Figure 14 Single Nozzle 60° Aeronhcll Model with Blunt Flow Interaction, 


Jarvinen and Adams Single Nozzle Results 




Non-dimensional Radial Coordinate 


• HLLE++, SSOR, and SST without compressibility correction 

• Literature says this is a steady flow, OVERFLOW results are unsteady 

• Overpredict shock standoff distances by ~ 20 % (used average distance) 

• Pressure comparison good for no jet case 

• For Ct =7, pressure is under-predicted at first data point and over-predicted at 

shoulder ^ 



Grid and Numerical Method Sensitivity 
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• HLLE++ with SST can predict the proper plume structure on this coarse 
Cartesian mesh 

• HLLC with SST and HLLE++ with SA incorrectly predict plume structure 
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Grid and Numerical Method Sensitivity (cont.) 
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• With a finer curvilinear mesh, HLLC with SST correctly predicts plume structure 

• HLLE++ with Spalart-Allmaras still predicts a steady “candle flame” behavior. 

• SST solutions unsteady 

* Shown above are instantaneous solutions in time and probably not at the same time 

• Running with SST compressibility correction or laminar makes solution more unsteady 
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Daso et al Single Nozzle Case 


NASA Marshall Trisonic 
Blowdown Wind Tunnel, 2007 

Apollo capsule with sonic nozzle 
Good Schlieren images 
No pressure data reported 
Freestream Conditions 
Mach=3.48, Re=4.88E6/ft 
Low thrust coefficient of 0.4 


A 
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Daso Results 


Daso et al (AIAA 2007-1423) 
Sonic nozzle 


■ 


Daso et al (AIAA 2007-1423) 
Sonic nozzle 


M = 3.48, C T = 0.4 Mach 0.0 0.5 1.0 1.5 20 2.5 3.0 3.5 4.0 4 5 5.0 


M = 3.48, C T = 0.4 


Density Based Shadowgraph 



I I 


Daso Results 


Daso 0.5" sonic nozzle 
Mass flow rate = 0.5 Ibm/s 
OVERFLOW Mach contours 
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• HLLE++, SSOR, and SST without compressibility correction 

• Great shock standoff distance comparison 

• Solution reaches a steady state 

• Makes this case good for grid refinement and numerical method 
sensitivity studies 12 


Grid Refinement Study 


• Axisym metric 

• Created script to generate mesh based on a global scaling 


parameter 

• Mesh spacing = value * scale factor 

1 The smaller the scale factor, the smaller the grid spacing, which 
means the greater the refinement 

Created 5 refinement levels 


Refinement 

Level 

Number of 
Points 

Increase per 
Direction 

1.25 

67581 

1 

1 

93051 

1.17 

0.5 

237867 

1.88 

0.25 

739590 

3.31 

0.15 

1767999 

5.11 i 



ref 0. 1 5 
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Qualitative Comparison to Schlieren 


• Features sharpen as 
refinement increases 

• Shock distance comparisons 
are more accurate with 
increased refinement 

• Feature position does not 
change much with last three 
refinements 

• Qualitatively, refinement 
level 0.5 is good enough 




ref 0.25 
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ref 0. 1 5 



Quantitative Comparisons 


ref 1 .25 



• Surface pressure coefficient on the capsule frontal area 

• Large gap between refinement levels 1 .0 and 0.5 

• Gap decreases with finer levels 

• Locations of local min/max appear to be converging 

• 3% change in drag coefficient between refinement levels 0.25 and 0.1 5 

• Negligible change in thrust coefficient after refinement level 0.5 

• Quantitatively, refinement level 0. 1 5 is best, 0.25 is probably good enough 
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Mach Along Centerline 


Daso Numerical Method Sensitivity (ref 0.25) 
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• All SST results (not using compressibility correction) are very similar 

• Correctly predict locations of the terminal shock, interface, and bow shock 

• Spalart-Allmaras results differ from SST results 

• Interface standoff distance 

• Pressure coefficient in recirculation area and near shoulder 

• SST with compressibility correction and laminar cases are both unsteady 

• SST without compressibility correction best choice 
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Jarvinen and Adams Triple Nozzle 
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Figure 49 Three Engine 40° AeroeheU Flow Field M * 2.0. C_ * 1.02. 0 * 10*. 

t TOTAL 


Three radially aligned nozzles 1 20 
degrees apart 

Modeled geometry behind the aeroshell 
as a solid piece 

Thrust coefficients of 0, 1 , 4, and 7 
Freestream Mach number of 2.0 



Jarvinen and Adams Triple Nozzle, Mach 2.0, Ct 7 


Problem Setup 



X (diameters) 


Used a curvilinear refined mesh for plume region 
Plenum Boundary Condition 

• BC4 1 - specify total pressure and temperature was unstable 
for this configuration 

• Adopted use of BC43 and BC3 1 - prescribed Q variables 
coupled with characteristic condition based on Reimann 
invariants 

• Started subsonic plenum region at the converging section of 
nozzle to encourage acceleration towards the nozzle exit 
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Jarvinen and Adams Triple Nozzle Results 



Model Front Face 


4)=0° 



Reasonable agreement for Ct= I 

Not very good agreement for Ct=4 

Okay agreement except at the nose for Ct=7 

For all three thrust coefficients, code-to-code 
comparison was much better than code-to-test 
comparison 
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Langley UPWT Pre-Test 


• Tunnel test designed for CFD 
validation 

• Used CFD to predict effects of 
model diameter 

• Wall effects 

• Possible blockage 

• Tunnel un-starts 

• CFD Run conditions 

• Mach numbers of 2.4, 3.5, and 4.6 

• Thrust coefficient range up to 10 

• Angle of attack sweeps up to 1 0 
degrees 

• 0, 1 , 3, and 4 nozzles 

• Experiment completed July 3 1 , 20 1 0 
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Pre-Test Results on Model Diameter 


• Modeling Assumptions 

• No attach hardware 

• Inviscid tunnel walls 

• Symmetric test section 

• Ran simulations of 4 and 6 inch model diameters 

• Saw notable wall effects for 6 inch diameter 

• 4 inch diameter too small for instrumentation 

• Selected 5 inch model diameter 


FUN3D (Bil Kleb, NASA LaRC) simulations for 6-inch diameter for C T =I0. 


-p. piii -I 0.6 076 0.1 0.06 0.01 0.01 046 at 076 0.6 I 



M=4.6 M=3.5 M=2.4 
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5 inch Pre-test Results 



WT pretest. 5 in diameter. M =2.4. 4 nozzle, C.= 10 


WT pretest. 5 in diameter. M =2.4. 1 nozzle. C,=5 
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Grid 

One focused plume grid for all nozzle configurations 

Overlap between plume and background grids cause blip in 
the bow shock 

• Ran a refined background case and saw minimal change in 
bow shock reflection location and model surface pressure 

Solutions 

Ran a series of cases representative of the test run matrix 
Results used for real-time test-to-code comparisons 
Preliminary comparisons to tunnel test results are promising 
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Pre-test Code-to-Code Comparison 
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•Good agreement for this case 
•Not all cases agree this well 

•Overset DPLR. used for this case 
(see next slide) 


23 


Overset DPLR 






diameters 


This work is from Kerry Trumble (NASA ARC) 

Point-matched grids are limiting for these geometries 

Made overset grid with Gridgen 

Domain connectivity with Suggar 

Will be using Usurp for force and moment calculations 
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Current Work: Post-Test Analysis 



• The recent Langley UPWT test provided a lot of well defined data 

• Current work for CFD team 

• Match the tunnel data 

• Explore effects of tunnel artifacts such as attach hardware and viscous tunnel walls 
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Future Work 
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TUNNEL STATION, in. 


ALL DIMESStOXS ARE /N INCHES 


9x7 Supersonic Wind Tunnel Modd Installation 


• OVERFLOW future work 

• More grid sensitivity study 

• Refinement, alignment 

• OVERFLOW grid adaption capabilities 

• Effects of turbulence modeling 

• Thermally vs. calorically perfect simulations 

• SRP team future work 

• Use the same model from the Langley UPWT 4’x4’ in the Ames 9*x7* tunnel 

• Hopefully obtain higher thrust coefficients 

• Less tunnel artifacts 


26 







Acknowledgements 


ETDP EDL Project 
SRP CFD and Wind Tunnel Teams 
JSC Applied Aeroscience and CFD Branch 


References 


• Jarvinen, P.O. and Adams, R.H., “The Aerodynamic Characteristics of Large Angled Cones with 

Retrorockets,” NASA Contract No, NAS7-576, Cambridge, MA, Feb. 1970 

• Daso, E. et al, “Dynamics of Shock Dispersion and Interactions in Supersonic Freestreams with 
Counterflowing Jets,” AIAA Journal Vol. 47, No. 6. June, 2009. 

• Trumble, K.A., Schauerhamer D.G., Kleb.W.B., Carlson, J.R., Buning, P.G., Edquist, K.T., Barnhardt, 

M.D.,“An Initial Assessment of Navier-Stokes Codes Applied to Supersonic Retro-Propulsion, 
AIAA 20 1 0-5047, June, 20 1 0. 


